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ABSTRACT 

The attenuation of starlight by interstellar dust is investigated in a sample of low redshift, disk-dominated 
star-forming galaxies using photometry from GALEX and SDSS. By considering broadband colors as a func- 
tion of galaxy inclination we are able to confidently separate trends arising from increasing dust opacity from 
possible differences in stellar populations, since stellar populations do not correlate with inclination. We 
are thus able to make firm statements regarding the wavelength-dependent attenuation of starlight for disk- 
dominated galaxies as a function of gas-phase metallicity and stellar mass. All commonly employed dust 
attenuation curves (such as the Calzetti curve for starbursts, or a power-law curve) provide poor fits to the ul- 
traviolet colors for moderately and highly inclined galaxies. This conclusion rests on the fact that the average 
FUV-NUV color varies little from face-on to edge-on galaxies, while other colors such as NUV-m and u — r 
vary strongly with inclination. After considering a number of model variations, we are led to speculate that 
the presence of the strong dust extinction feature at 2 175 A seen in the Milky Way (MW) extinction curve is 
responsible for the observed trends. If the 2175A feature is responsible, these results would constitute the first 
detection of the feature in the attenuation curves of galaxies at low redshift. Independent of our interpretation, 
these results imply that the modeling of dust attenuation in the ultraviolet is significantly more complicated 
than traditionally assumed. These results also imply a very weak dependence of the FUV-NUV color on total 
FUV attenuation, and we conclude from this that it is extremely difficult to use only the observed UV spectral 
slope to infer the total UV dust attenuation, as is commonly done. We propose several simple tests that might 
finally identify the grain population responsible for the 2175A feature. 
Subject headings: dust, extinction — ultraviolet: galaxies — galaxies: ISM 



1. INTRODUCTION 

Dust is everywhere. Its obscuring and emissive signa- 
tures have been observed in environs ranging from nearby 
dwarf irregulars to giant ellipticals, within our own and 
nearby galaxies, in intergala ctic space, and even in the 
highe s t redshift quasars ( e .g.. lTrumplerlll930l iDisnev et alJ 
1989 t iKnapp et all IT989t iGoudfrooii et alj 11994 lZaritsk\| 
1994 1 IXu & Buatlll995tlWang & Heckman|ll996tlElbaz et all 
19991 IXilouris et al.lll999t ICalzettill2001h Motta et alj|2002t 
Wang et alj|2008t iMenard et alj|2009h . 



Our understanding of the composition and spatial distri- 
bution of dust in a variety of environments has increased 
markedly in the past decade, thanks largely to the launches 
of the Spitzer Space Telescope and the Galaxy Evolution Ex- 
plorer (GALEX). For example, the fraction of a galaxy's gas 
mass that is locked up in dust is now well-characterized as a 
function of galaxy mass, star fo r mation rate (SFR), and metal - 
licity (e.g.. iDraine etall 120071: Ida Cunha et al]|2008l l20Tob . 
The radial distributi on of gas and dust has also been ex- 
tensively studied (Munoz-Ma teos et al.l 120091) . The impor- 
tance and prevalence of polycyclic aromatic hydrocarbons 
(PAHs) is now widely acknowledged tha nks to moderate reso- 
lution spectroscopy in the mid-IR (e.g . , iLeger & Puget|[T984l : 
iSmith et alJl2007tlO'Dowd eTaTl2009l) . 

Despite these advances, basic questions remain. It is still 
not known, for example, if the majority of dust forms in AGB 
outflows, in supernovae ejecta, or condenses directly out of 
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the ISM (Gallia no et a"T]|2008l:lDrainell2009l) . Computation of 
the dust extinction curve from first principles relies on a num- 
ber of uncertain inputs, including the size, shape, and com- 
position of the grains, and the optical constants of the con- 
stituent materials (e .g., Weingartner & Draine 2001; Draine 
120031: iGordon et al]|2003l IZubko et all 120041 ) A further un- 
certainty with broad implications for our understanding of 
galaxies is the relation between the stellar and dust content of 
galaxies; for example, how the grain size distribution varies 
within and between galaxies. Due to gas inflow and outflow, 
the stellar and gas-phase metallicities are not simply related, 
further compounding the problem. Without a solid theoreti- 
cal understanding of the complex relation between stellar and 
dust content, these two components need to be modeled flex- 
ibly and independently. Despite this widely acknowledged 
fact, in practice a particular dust obscuration model is often 
simply assumed in order to then infer the underlying stellar 
population. Such an assumption can introduce significant bi- 
ases, as we will discuss below. 

The net effect of dust obscuration within a galaxy is the 
result of the combination of the underlying dust extinc- 
tion curve, geometry, and radiative transfer. Here, geom- 
etry is taken to mean not only the large-scale distribution 
of dust with respect to the stars, but also the local geome- 
try, which includes the dumpiness of the interstellar medium 
(ISM). Geometrical effects can make the relation between 
the underlying extinction curve and the resulting attenua- 
tion curve (which measures the net l oss of photons within a 
galaxy) arbitrari l y complicated (e.g., iNatta & Panagial 1 19841 ; 
[ Witt et all [1991 ICalzetti et al.l Il994t IWitt & Gordod Il996t 
Varos i & Dwekl 119991 IWitt & Gordon! I2000t iGranato et all 
20001: ICharlot & Falll l2000t iPierini et alj 12004 iTuffs et all 
12004 ; lPanuzzo et alJ2007l) . For example, a clumpy ISM, aris- 
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ing for example from turbulence, will result in an attenuation 
curve that is greyer than the extinction curve. 

When the optical depth approaches unity, the effects of ge- 
ometry become especially important. In such cases, the light 
measured in the blue/ultraviolet may not trace the same re- 
gions of the galaxy as the red/near-IR light. Consider a disk 
galaxy viewed edge-on that is optically thick in the UV and 
optically thin in the near-IR. The ultraviolet light emitted from 
the far side of the galaxy will be heavily extinguished, and so 
the ultraviolet light collected by an observer will be prefer- 
entially sampling the near side of the galaxy. In contrast, the 
near-IR will more faithfully trace the entire stellar population 
because the galaxy is relatively transparent in this wavelength 
range. These effects are rarely considered when modeling the 
SEDs of galaxies. 

If we knew the underlying, unobscured stellar population, 
the ratio between the intrinsic and observed light would then 
provide a direct measure of the attenuation. Unfortunately, it 
is not possible to know unambiguously the underlying stel- 
lar population. Less direct methods are therefore required. 
For example, the attenuation curve can be estimated in a sam- 
ple of galaxies with similar stellar populations and variable 
amounts of dust. In such cases, ratios between more and less 
heavily attenuated spectra can provide estimates of the net at- 
te nuation curve in the s ample. Such a technique was utilized 
bv lCalzetti et al.l d 19941) in order to probe the average attenu- 
ation curve of starburst galaxies. The Balmer decrement was 
used as an estimator for the amount of dust attenuation, and, 
with the assumption that the galaxies in their sample were of 
similar metallicity and SFR, the average attenuation proper- 
ties were estimated. 

Similar techniques for probing the wavelen gth-dependent 
attenua tion have recently been applied by [Johnson et alj 
(l2007allbh to photometr y of a large sample of low redshift 
galaxies. lConroyl (120 1 01) who analyzed the restfram e UV pho- 
tometr y of star-forming galaxies at z ~ 1, and iNoll et aD 
(2009) who analyzed spectra of star-forming galaxies at z ~ 2. 
In the present work we follow in this vein by considering 
the colors of disk-dominated galaxies as a function of their 
inclination. Since inclination will correlate with dust atten- 
uation but not with stellar populations, comparing less to 
more inclined systems allows us to isolate the effects of in- 
creasing dust opacity on the observed properties of galax- 
ies (e.g., iGiovanellietail 11994 [T995; Masters et alj 2003; 
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Another constraint on the attenuation curve comes from 
the relation between the ratio of total infrared to UV lumi- 
nosity and UV spectral slope (Ljir/^uv _ P, or the TRX- 
/?' relation). Star-forming galaxies with redder UV spec- 
tra tend to have higher IRX. It is widely believed that this 
relation is prim arily a sequence in dust attenuation (e.g., 
iKong et aTll2004l) . Various recipes have thus been proposed 
to use the IRX-/3 relation to estimate the total U V attenua- 
tion based on the observed UV spectral slope (e.g . , iBuat et ail 
l2005t iBurgarella et al.ll200l ICortese et al.ll2008l) . This rela- 
tion is of particular importance to the study of high-redshift 
galaxies, where only restframe UV and optical photons can 
be readily collected for large samples of galaxies (although 
significant samples at z ~ 1-2 with res t frame IR data are 
growing rapidly, e.g., iReddv et al.1 120061: ISalim et al.l 12009b 
Redd v etalj 120101) . For most studies of high redshift galax- 
ies, the observed UV slope, f3, is used in conjunction with a 
locally estimated IRX-/3 relation to estimate the dust opacity. 



This approach is essential, for example , to interpret recent ob - 
servations of galaxies at z ~ 6-8 (e.g.. Bou wens et aT]| 2009). 

The IRX-/3 relation is different for starbu rsts dMeurer et al.l 
1999) and normal star-forming galaxies dDale et al.1 120071 : 
Boissier et aTl l2007h . and i n addition depend s somewhat on 
the sta r formation history (IKong et al.l l2004t iJohnson et alJ 
2007b : iMunoz-Mateos et al l 120091) 7 60/xra luminosity 
dTakeuchi et al.l 120101) . bolometric lu minosity (IReddv et al.l 
120061) . and the star-dust geometry dPanuzzo et al.l 120071) . 
These dependencies result in substantial scatter in the IRX-/3 
plane, which, in conjunction with the fact that the IRX-/3 
is nearly vertical over much of the relevant parameter space 
(i.e., IRX varies considerably over a relative ly narrow range 
in /3), calls its utility into question (Bell 2002). 

The most prominent feature of the MW extinction curve 
is the strong, broad dust feature at 2175A, the 'UV bump' 
dStecherll!965l) . This feature is also seen, albeit more weakly, 
along most sightlines in the L MC and along one of the five 
sightlines probed in the SMC dGordon et al.ll2003l). It is also 
seen along sightlines in M31 dBianchi et alJ 119961) . There 
are now confident detections of this feature in the extinc- 
tion curves of more distant galaxies as well, as probed by 
backgr ound gamma ray bursts and quasars, and gravitational 
lense s (Elfasdot tFet al J 120091: iMotta et al.1 120021: IWang et al.1 
l2004t [Mediavill a et al.ll2005l) . Puzzlingly, there are many ex- 
amples of galaxies t hat do not show this feature in the ir ex- 
tinction curve (e.g.. lYork et~aLll2006t IStratta et alJl2007l) . The 
carrier of this extinction feature is not known, although ow- 
ing to its strength it must be due to some abundant ma- 
terial, such as carbon. Measurements of the grain albedo 
suggest that th e feature is not due to scattering (see data 
compilations in lWitt & Gordonll2000tlDraindl2003h . Several 
dust models associate this fe ature with PAH absorption (e.g., 
Weing artner & Draind 120011) alth ough there are other possi- 
bilities dDraine & Malhot raTl993L e.g.,). 

Understanding the prevalence of the UV bump in the at- 
tenuation curves of other galaxies is essential for broadband 
photometric studies of galaxies. For example, at z ~ the 
UV bump falls into the GALEX NUV-band, at z - 1 it falls 
within the B-band, and at z ~ 2 it redshifts into the /?-band. 
If present in the attenuation curves of galaxies, the UV bump 
would therefore result in substantially more, and likely more 
uncertain, attenuation in particular filters at particular red- 
shifts. 

Based on spectra from the International Ultraviolet Ex- 
plorer, Calzetti et al. found an average attenuation curve 
for starburst galaxies that lacked a UV bump. The ab- 
sence of a UV bum p in the starburst attenuation curve led 
IWitt & Gordonl d2000l) to suggest that in such galaxies the un- 
derlying dust extinction curve lacks a UV bump. A detailed 
analysis of the ultraviolet through infrared photometry of M5 1 
revealed little evidence for a UV bump within individual HII 
regions dCalzetti et all 120051) . IConrovl d2010h also found no 
evidence for a UV bump as strong as that seen in t he MW in 
star-fo rming galaxies at z ~ 1 . In stark contrast, INoll et all 
(2009) presented strong evidence for a UV bump in stacked 
restframe UV spec tra of z ~ 2 star-forming galaxies. Finally, 
Capak etafl (120091) presented tantalizing evidence for a strong 
UV bump in at least one z ~ 7 galaxy. 

Thus, while the UV bump appears to be a ubiquitous feature 
of the MW and LMC extinction curves, there is little evidence 
for this feature in the net attenuation curves of other galaxies. 
However, it is important to recognize that a systematic inves- 
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tigation of the UV attenuation in z ~ 'normal' star-forming 
galaxies is currently lacking. While the effects of geome- 
try and scattering can diminish the strength of the UV bump 
with respect to the underlying extinction curve, it is a generic 
prediction of radiative transfer calculations that if the UV 
bump is present in the extinction curve, its presence will be 
detectable in the resulting attenuation curve dWitt & Gordon! 
l2000tlPanuzzo et al.ll2007l) . Given the observations of the UV 
bump in the extinction curves of other galaxies, the attenua- 
tion curves of at least some galaxies should show a UV bump, 
even if the strength of the bump is weak. This expectation 
serves as motivation for the present study. 

There is currently no clear picture linking together these 
various observations. One possibility is that the radiation 
fiel d modulates the stre ngth of the UV bump, as suggested 
by iGordon et al.1 (120031) . The typical radiation field in the 
SMC and the local starburst galaxy sample is much harsher 
than the MW, potentially explaining why a UV bump is not 
seen in such systems. It would be hard to explain the results 
of Noll et al. in this context, however, given that the i r sam- 
ple is dominated by high SFR systems. iGranato et al.l (2000) 
explained the result from Calzetti et al. as being due to the 
fact that the UV energy production in starbursts is dominated 
by young stars that are heavily embedded within molecular 
clouds. In their model, attenuation in starbursts is therefore 
governed by the wavelength-dependent fraction of UV pho- 
tons produced by young stars (whose light is heavily extin- 
guished by their birth cloud). These authors predicted that the 
attenuation suffered by 'normal' star-forming galaxies would 
show evidence of a UV bump because in this case significant 
UV energy is provided by intermediate age stars that have 
left their birth clouds. Such stars will suffer attenuation pri- 
mary from the diffuse dust, where significant 2 175 A absorp- 
tion may be expected. It is not immediately obvious whether 
or not this explanation can ac commodate the results from 
Noll et al. iDraine et al.l (120071) has recently demonstrated a 
deficiency of PAH emission in the i nfrared in galaxies with 
low gas-phase m etallicities (see also Engelbracht et al.ll2005t 
Smith et al. 2007). Metallicity may therefore play an impor- 
tant role. 

In the present work we investigate the wavelength- 
dependent attenuation by dust for a sample of disk-dominated 
galaxies. Our sample is carefully selected to be homoge- 
neous and complete. We then consider the ultraviolet, op- 
tical, and near-infrared colors as a function of inclination. 
Since inclination will only correlate with dust attenuation and 
not physical parameters such as star formation nor metallic- 
ity, the inclination-dependent colors will provide a robust and 
sensitive probe of the wavelength-dependent attenuation in 
disk-dominated s ystems. This technique is t herefore similar 
in spirit to that of iCalzetti et all (11994 12000), although for a 
sample of 'normal' star-forming galaxies, and utilizing pho- 
tometry rather than spectroscopy. Averaging the flux of many 
galaxies considerably simplifies the modeling of the underly- 
ing stellar population because the average star formation his- 
tory (SFH) of many normal galaxies must be smooth. By nor- 
mal here we mean galaxies not chosen to have special SFHs, 
such as starburst or post-starburst galaxies. We will focus es- 
pecially on the attenuation properties in the NUV band, and 
will therefore be able to make strong statements regarding the 
presence of the UV bump in the average attenuation curve of 
normal star-forming galaxies. 

Where necessary, a flat ACDM cosmology with (fi m , Oa) = 



(0.30,0.70) is adopted, along with a Hubble constant of Hq = 
100/i km s" 1 Mpc" 1 . A ll magnitudes ar e in the AB system 
dOke & Gunnlll983l) . A lChabrierf d2003l) initial stellar mass 
function (IMF) is adopted when quoting stellar masses. 

2. DATA 

Redshifts, coordinates, ugriz photometry, and structural 
properties are de rived from Sloan Digital Sky Survey (SDSS; 
lYorket al.1 120001) data, as made availa ble through the NYU 
Value Added Galaxy Catalog^ (VAGC; IBlanton et all 120051) . 
We do not simply use the photometry from the SDSS cat- 
alogs, because the automated pipeline often incorrectly de- 
blends large galaxies. In addition, the images themselves are 
generally over-subtracted due to the sky estimation procedure, 
which uses a median smoothed field in 100" by 100" boxes 
( IBlanton et al.ll2005l) . In order to address these issues, Blan- 
ton et al. (in prep) measure the smooth sky background from 
heavily masked versions of each imaging run. In this way, we 
avoid subtracting away light from large galaxies. Blanton et 
al. (in prep) have im plemented a de-blen ding pipeline similar 
to that described by iLupton et al.l (1200 ll) . but with parame- 
ters optimized for brighter and larger objects. For small ob- 
jects these measurements agree well in general with the orig- 
inal SDSS catalog values; for larger objects (half-light radii 
> 10") our galaxy measurements tend to be brighter in flux 
and larger in size. 

We use the 'lowz' VAGC sample, which is constructed from 
the SDSS Data Release 6 dAdelman-McCarthv et al] 120081) . 
but with the photometric re-reductions discussed abov e. This 
sample includes all galaxies with redshifts z < 0.05. ISersid 
( 1968) profile fits have been made to the radial light profiles 
of galaxies in this sample (the profile fits were corrected for 
seeing). We will make use of the Sersic index n, where recall 
that n = 1 corresponds to an exponential light profile and n = 4 
corresponds to a de Vaucouleurs profile. We will also make 
use of the axial ratio, b/a, where b and a are the semi-minor 
and semi-major axes, respectively. The axial ratio is measured 
within a radius containing 90% of the total light, and will be 
taken to be a proxy for the galaxy inclination. Petrosian mag- 
nitudes are used herein. 

Gas-phase oxygen abundances have been estimated from 
the SDSS spectra by iTremonti et al.1 ((2004), although these 
measurements are only available for galaxies in the SDSS 
Data Release 4. When considering galaxies split according 
to gas-phase metallicities, we will only use galaxies that are 
in the Tremonti et al. catalog. 

Near-IR photomet ry are provided by the Two Micron All 
Sky Survey (2MASS lJarrett et al.ll200ot) . We make use of the 
/T-band only in Figure Q] when investigating the properties 
of our sample, and so make no attempt to perform matched 
aperture photometry on the 2MASS data. 

The GALEX survey dMartin et al.ll2005l) is providing ultra- 
violet imaging in the FUV and NUV filters (with effective 
wavelengths of 1539A and 2316A, respectively). In order to 
accurately probe the UV through optical SED of galaxies as 
a function of inclination, we have re-measured GALEX pho- 
tometry self-consistently for our sample by computing magni- 
tudes within the Petrosian aperture defined by the SDSS pho- 
tometry (Petrosian magnitudes in SDSS are defined accord- 
ing to the Petrosian radius, rp, measured in the r-band). Our 
photometry was performed on GALEX images that have been 
calibrated and background subtracted. Thus, the FUV, NUV, 

4 http://sdss. physics . nyu . edu/vagc/ 
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FIG. 1. — Properties of the low redshift galaxy sample used in the main 
analysis. Upper left: Redshift distribution normalized to the total number 
of objects in the sample. Upper right: Relation between Sersic index n and 
axial ratio b/a for the full low redshift sample (randomly diluted by 60% 
for clarity; grey symbols), compared to the subsample used in the main anal- 
ysis (black symbols). Lower left: Histogram of axis ratios for the sample 
used in the main analysis. Bin widths are 0.1, starting at b/a = 0.2, which 
corresponds exactly to the bins used in later sections. The hatched bin is 
excluded from our analysis for reasons discussed in the text. Lower right: 
Color-color diagram (symbols are as in the upper right panel). UV-optical- 
near-IR colors have been shown to efficiently separate truly quiescent from 
dusty star-forming galaxies. The galaxies used in this sample populate the 
latter locus. 

and ugriz magnitudes are all measured through the exact same 
aperture. The Petrosian aperture (which is defined as 2rp) is 
much larger than the half-light radius and so differences in the 
PSF between GALEX and SDSS should have a negligible ef- 
fect on the derived photometry, especially for the exponential 
light profiles that dominate our sample. Ze ro points for the 
FUV a nd NUV magnitudes are adopted from Morriss ev et alj 
d2007h . 

E(B - V) reddening estimates are derived from 
Schlegel efaD (fl998h via the utilities in the kcorrect 
v4 . 1 . 4 software package (Blanton & Roweis 2007). Pho- 
tometry is corrected for Galactic extinction using E(B — V) 
and the MW extinction curve of lCardellietaT] ([1989) . In 
particular, we use the following values for = Aj/E(B — V), 
where A, is the extinction in magnitudes for filter z: 8.29, 
8.18, 5.16, 3.79, 2.75, 2.09, 1.48, and 0.30 for the FUV, NUV, 
ugriz, and K bands. 

K-corrections are estimated with the routines made avail- 
able in the kcorrect v4 . 1 . 4 package. The SDSS ugriz 
photometry is used as input to estimate k-corrections for both 
the SDSS photometry and GALEX photometry. Photometry 
is k-corrected to z = 0.04, which is the median redshift of the 
sample. Typical k-corrections are of order 0.01 magnitudes 
and therefore have no significant effect on our results; they 
are included merely for completeness. 

The kcorrect v4 . 1 . 4 package also provides estimates 
of the stellar masses of galaxies given a set of input photom- 



etry. We have used the SDSS photom etry to es t imate stel- 
lar masses for our sample, assuming a Chabrieri (120031) ini- 
tial mass function. We prefer to use photometry for stellar 
mass estimation because masses derived in this way are only 
weakly dependent on galaxy inclination (stellar masses differ 
by < 0.1 dex betw een highly inclined and face-on systems; 
iMaller etail 120091) . O ther, popular stellar mass estimators 
based on spectroscopy ([Kauffmann et al. 2003; Toieir o et"aT] 
2009), suffer stronger dependencies with inclination, which 
indicates a systematic source of bias in the masses. 

Our goal is to explore the inclination-dependent dust char- 
acteristics of a uniform sample of disk-dominated galaxies. 
In order to construct such a sample we select galaxies with 
n < 2.5 to separate disk-dominated from bulge-dominated 
galaxies. In addition, galaxies are selected in a narrow stel- 
lar mass range 9.5 < log(M/M Q ) < 10.0 in order to de- 
fine a homogeneous and volume-li mited sample. Mass and 
SFR are correlate d at z ~ (e.g., Brinch mann et al.l I2004t 
ISalim et ail 12007b and so our selection on mass may be in- 
terpreted roughly as a selection on SFR as well. Galaxies at 
z < 0.01 are removed from our sample, as they are generally 
very large on the sky and thus accurate photometry is compli- 
cated. We also remove 112 galaxies with rp < 4", as inclina- 
tion measurements become unreliable for such small galaxies 
because of seeing (see below). With these cuts we are left 
with 3394 galaxies with detections in all GALEX and SDSS 
filters. This is our fiducial sampl e. 

As discussed in iMasters et al.l (120031) . the effects of seeing 
can bias inclination measurements for galaxies with small an- 
gular sizes. The bias manifests itself as a deficiency of small, 
highly inclined objects. Indeed, our sample shows such a de- 
ficiency for galaxies with rp < 8". The median Petrosian ra- 
dius of our fiducial sample is 8". 4, so this bias potentially 
affects a significant fraction of our sample with high incli- 
nations. However, we have re-computed all of the results 
in following sections only including galaxies with rp > 8" 
(and widening the stellar mass cut to 9.25 < log(M/M Q ) < 
10.0 in order to increase the sample size), and find that the 
inclination-dependent colors change by an inconsequential 
amount. Moreover, these issues primarily affect galaxies with 
0.2 < b/a < 0.3 and since we consider only galaxies with 
b/a > 0.3 we can be confident that our results are robust to 
seeing biases. 

Several properties of our fiducial sample are shown in Fig- 
ure Q] including the redshift distribution, scatter plot of in- 
clination and Sersic index, distribution of inclinations, and a 
color-color plot. The deficiency of galaxies with b/a < 0.3 
is due in large part to the seeing bias noted in the previous 
paragraph. However, one would also expect the intrinsic num- 
ber of axial ratios to drop at such low values because of the 
rareness of intrinsically thin disks. This bin is excluded from 
our analysis below. 

The NUV-r color, in combination with an optic al-near-IR 
color such as r-K, cleanly separates star-forming and quies- 
cent galaxies, even when t here are substantial amounts of dust 
in sta r-forming galaxies (Willia ms et al J 120091 iBundv et all 
2009). It is clear from this figure that our sample is domi- 
nated by star-forming galaxies, which is not surprising given 
the mass range of the sample. Galaxies in this mass range 
have current SFRs that are comparable to their past-averaged 
SFRs (i.e., the birth-p arameter for these objects is ss 1; 
iBrinchmann et a l. 2004). In other words, star formation has 
proceeded in an approximately continuous manner over the 
lifetime of these galaxies. The results presented in following 
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sections are unchanged if we remove the few truly quiescent 
galaxies in our sample (galaxies with NUV-r colors greater 
than « 4.5). 

2.1. On the challenges and benefits of averaging many 
galaxies 

In Sj4] we will present average colors of disk-dominated 
galaxies as a function of inclination. Within each bin in in- 
clination the average colors are computed by averaging the 
fluxes for typically several hundred galaxies. In the absence 
of dust, averaging in this manner is equivalent to construct- 
ing a 'super-galaxy' that is composed of all of the stars of the 
galaxies in the bin. With the addition of dust attenuation, the 
interpretation of this averaging procedure is more complex. 
In this case magnitudes may be somewhat simpler to inter- 
pret because magnitudes are linear in the attenuation. Indeed, 
if the underlying stellar population is unchanged and only the 
dust content varies, then averaging magnitudes is probably the 
correct approach. In ij4]we will present average colors com- 
puted both by averaging fluxes and averaging magnitudes in 
order to demonstrate that our results are insensitive to this dis- 
tinction. Unless explicitly stated otherwise, our results will be 
computed by averaging fluxes. 

The major benefit of averaging many galaxies is that the 
SFH of the ensemble must be smooth. As noted in the intro- 
duction, this need not necessarily be the case if one selects 
special classes of galaxies such as starburst or post-starburst 
galaxies. Our sample is composed of normal star-forming 
galaxies, and so our assertion that the SFH must be smooth 
is justified. One thus cannot appeal to sharp variations in the 
SFH (e.g., recent bursts) to explain the average colors of the 
stacked galaxies. This will prove to be an important asset 
when interpreting the inclination-dependent colors in £|4] 

3. MODELS 
3.1. Stellar population synthesis 

The stellar population sy nthesi s (SPS ) treatment closely 
follows that o flConroy etail (120091 l20Toh . to which the reader 
is referred for details. Briefly, the code produces the evolu- 
tion in time of the spectral energy distribution of a coeval set 
of stars from 10 6 5 to 10 10 15 yr, for a wide range in metal- 
licities. Stellar evolution from the main sequence through the 
thermally- pulsating AGB an d post-AGB phases are included. 
The I MF oflChabrier|(|2003l) i s adopted. 

In IConroy & Gunnl (1201 Oh . this SPS code was exten- 
sively calibrated against a suite of observations and was also 
compared to several other popular SPS codes. The ver- 
sion of t he SPS code use d here in is identical to that dis- 
cussed in IConroy & Gunnl (120101) . We have also compared 
this SPS code to the predictions of the Starbu rst99 code 
dLeifherer et al.l[l999t IVazquez & Leifhererl l2005) for a con- 
stant SFH at solar metallicity. The spectra from these two 
models are very similar in the ultraviolet and optical. This 
agreement is encouraging because Starburst99 contains a 
more sophisticated treatment of massive star evolution and 
also includes nebular continuum emission, which is not in- 
cluded in our models. 

The underlying stellar population in our default model will 
consist of solar metallicity stars that formed in a continuous, 
constant mode of star formation from t = 0.0 to t = 13.7 Gyr. 
In order to explore the impact of these assumptions we will 
also consider models with metallicities Z = 0.5Z Q and Z = 
1.5Z Q , and models where a fraction, c < 1, of mass forms 
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FIG. 2. — UV attenuation curves, normalized to the attenuation in the 
V— band. The MW curve with Ry = 3.1 is repeated in all panels for refer- 
ence. Top left: The MW curve (with Ry = 3.1) both with and without a UV 
bump are compared to the Calzetti et al. attenuation curve and a MW curve 
with Ry = 2.0. The transmission curves for the FUV, NUV, and H-band filters 
are also shown. Top right: Power-law curves of the form t(A) oc A , with 
8 = —0.4, —0.7, —1.0, and —1.3. Bottom panels: Attenuation curves derived 
from the models of WG00 for MW and SMC extinction curves are shown 
both for the 'shell' and 'dusty' geometries. The V-band normalized atten- 
uation curve for the dusty geometry depends on the dust column density in 
the sense that larger columns yield grayer curves. In contrast, the shape of 
the attenuation curve for the shell geometry is independent of column density 
and so only one curve is seen for that geometry. 

in a constant mode of star formation, the rest having formed 
instantly at t = 0.0. All colors considered herein have been k- 
corrected to the median redshift of the data, which is z = 0.04. 



3.2. Dust models 

A suite of phenomenological dust models are considered 
for the attenuation of starlight. The models are similar in 
spirit to the two-component dust model of ICharlot & Falll 
(2000). This model distinguishes between attenuation suf- 
fered by young stars in their natal clouds and attenuation of 
all starlight due to diffuse cirrus dust. 

Attenuation due to diffuse cirrus dust is described by the 
V-band opacity Tyj = tj and an attenuation curve. By 'at- 
tenuation curve' we mean the variation in attenuation opti- 
cal depth with wavelength, normalized to 1.0 in the V-band, 
i.e., t(X)/tv- In this model, stars with ages younger that f esc 
are subject to additional attenuation (attributed to their birth 
cloud) characterized by an optical depth Ty.bc = i~b c with an 
attenuation curve that in principle may differ from the curve 
describing the diffuse dust. Based on data from a sampl e of 
low redshift star-forming galaxies, Ch ariot & Falll d2000l) fa- 
vor values of Tb c =1.0 and f esc = 10 7 yr with an attenuation 
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curve of the form r oc A -0 7 for both components. Herein we 
will consider additional attenuation curves and model param- 
eters, although for our default model we will adopt r\, c = 0.5 
and f esc = 10 7 yr. 

This two component model has strong observational mo- 
tivation not only from direct observations of young stars em- 
bedded in molecular clouds but also from integrated spectra of 
star-forming galaxies, where the opacity measured in balmer 
emission lines is a factor of roughly two larger than the opac - 
ity measured from the stellar continuum (ICalzetti et al.l 1994b . 

We will also briefly consider a dust model that, in ad- 
dition to the two component model described above, al- 
lows for a 'skin' of completely unobscured starlight, quan- 
tified as the fraction, /, of the stellar mass that is unob- 
scured. Motivation for such a model comes from recent ob- 
servations of extend ed UV emission in the outer edges of 
many disk galaxies (Thil ker et al.ll2007h . The young stars in 
these outer regions are apparently metal-poor (Z ~ Z Q /10), 
and the extinction they suffer is smal l, though non-zero 
dGil de Paz et al.ll2007t PWerk et al.ll2010b . In addition, it is 
well-known that the d ust-to-gas ratio decreases with radius 
in disk galaxies (e.g.. Ilssa et al.l [19901 iBoissier et al.l I2004t 
M unoz-Mateos et alj 2009), and so the outer regions of galax- 
ies should suffer substantially less extinction than the inner 
regions. Finally, a sizeable fraction (« 10-30%) of O B stars 
are 'r unaways' that have escaped their birth clouds (IStond 
1991). These stars will likely still suffer attenuation from the 
diffuse dust component, but nonetheless they may constitute 
another skin-like population that contributes a bluer than av- 
erage spectrum to the integrated light. 

3.2.1. Attenuation curves 

We now summarize the various attenuation curves that are 
considered herein. Simple power-law curves are considered, 
characterized by an index S such that r oc X s . We also con- 
sider attenuation curves whose functional for m is equivalent 
to the average extinction curve in the MW (Carde iliet al.1 
1989). As mentioned in the Introduction, attenuation and ex- 
tinction are different concepts, but in simple geometries they 
can be roughly equivalent (modulo wavelength dependence of 
the dust albedo) Q When considering the MW curves we allow 
for the freedom to vary the strength of the UV bump, param- 
eterized by B. In the Appendix we provide formulae for con- 
structing MW-like extinction curves with arbitrary UV bump 
strength. The default MW curve is recovered for B = 1.0; for 
B = 0.0 there is no bump. Below we will consider an attenu- 
ation curve that resembles the MW with B = 1 .0, B = 0.8, and 
B = 0.0. We also allow for the freedom to vary the standard 
extinction parameter Ry = Ay / 'E(B-V). Note that a smaller 
Ry corresponds to a steep extinction curve. The average ex- 
tinction curve in the MW is well-fit with Ry = 3.1, while a 
value of Ry w 2 provides a good fit to the average SMC ex- 
tinction curve at A < 5000A (|Pellll992l) . Below, when we 
refer to 'MW extinction' generically, we mean the standard, 
average MW curve with Ry = 3.1 and B = 1.0. 

Attenuation curves produced by the models of 
iWitt & Gord"onl ( 120001 WG00) will also be considered. 
These authors generate attenuation curves from both MW 
(adopting Ry = 3.1) and SMC extinction curves, taking 
into account the effects of scattering and geometry via 
radiative transfer calculations. They produce curves for three 

5 Later, when we refer to 'MW attenuation' we mean an attenuation curve 
that is equivalent to the observed, average MW extinction curve. 



large-scale star-dust geometries: a 'shell' geometry where 
all of the dust is in front of the stars (i.e., a uniform screen), a 
'dusty' geometry where the stars and dust are equally mixed 
within a sphere, and a 'cloudy' geometry where stars and dust 
are equally mixed within the inner 70% of the system, with 
the remaining outer portion containing stars but no dust. The 
WG00 models assume spherical symmetry. For each of these 
large-scale geometries WG00 consider both a homogeneous 
and clumpy distribution for the local distribution of dust. 
There are twelve models in all (two extinction curves, three 
large-scale geometries, and two local geometries). Each 
model is computed for a wide range in column densities. 

The attenuation curves considered herein are shown in Fig- 
ure |2] For the WG00 models, only the homogeneous local 
geometries are shown, and only for the shell and dusty large- 
scale geometries. Notice that the WG00 model constructed 
with a MW extinction curve and a shell geometry produces a 
UV bump that is stronger than is observed in the MW extinc- 
tion curve. This follows from the observationally-motivated 
assumption that the UV bump is a true absorption feature. In 
contrast, the dusty geometry produces bump strengths both 
stronger and weaker than the bump in the MW extinction 
curve, depending on the overall opacity. This serves to high- 
light the fact that UV bumps can be either stronger or weaker 
than the underlying extinction curve depending on the large- 
scale star-dust geometry. Notice also that the dusty geome- 
tries produce attenuation curves that are functions of the total 
dust column density in the sense that higher column densities 
produce grayer V-band normalized attenuation curves. 

In this work we will consider the variation in broadband 
colors of galaxies as a function of inclination. Based on our 
understanding of disk-dominated galaxies, we expect the pa- 
rameter tj, which characterizes the V-band opacity of the 
diffuse dust, to vary with inclination. The attenuation curve 
may also vary with inclination, becau se the dust column den- 
sity will vary with inclination (e.g.. i Varosi" &Dweki ri999: 
IWitt & GordoifeOOOt iTuffs et al.ll2004l) . However, we do not 
expect the dust attenuation suffered by young stars to depend 
on inclination, as this is a local phenomenon. We will there- 
fore not allow the parameter t\, c to vary with galaxy inclina- 
tion. For our main results we will also fix the attenuation 
curve for the young component to be a power-law with index 
S = -1 .0. Notice that there is no UV bump associated with the 
birth cloud dust. 

4. RESULTS 
4.1. UV colors 

In Figure [3] we compare the average UV colors of a sam- 
ple of disk-dominated galaxies as a function of inclination 
to a suite of models. Recall that the data sample consists of 
galaxies at 0.01 < z < 0.05 with Sersic index n < 2.5 and 
stellar masses in the range 9.5 < log(M/M Q ) < 10.0. Data 
are binned in inclination with width 0.1, starting at b/a = 0.3. 
Within each bin the average flux in each filter is computed. 
The resulting colors redden monotonically with increasing in- 
clination. Due to the large sample size, the errors on the mean 
colors are in all cases smaller than the symbol sizes. We have 
computed errors both via bootstrap re-sampling and the naive 
variance and find comparable results. 

The model predictions all have the same underlying stellar 
population, which consists of stars with solar metallicity hav- 
ing formed in a constant mode of star formation from t = 0.0 
to t = 13.7 Gyr. For each model starlight is attenuated ac- 
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FIG. 3. — UV color-color diagrams comparing stellar population models to the average colors of low redshift disk-dominated galaxies as a function of 
inclination. The data points become monotonically redder in NUV— u with decreasing b/a (increasing inclination). All models assume a constant SFH, are at 
solar metallicity, and adopt the same attenuation prescription for young stars. The models differ only in the adopted dust attenuation curve for the diffuse dust. 
For each attenuation curve except for those from WG00, the trajectory in color-color space is due to an increase in the diffuse dust opacity, tj, from 0.0 to 1.5. 
The opacity increases by intervals of 0.3, as indicated by the symbols along the models. Top left: Milky Way attenuation curves with Ry = 2.0 and with different 
treatments of the UV bump at 2175A: standard MW bump strength, no bump, 80% relative bump strength. Top right: Attenuation curves based on the WG00 
dust models. Models are constructed with both SMC (dashed lines) and MW (dot-dashed lines) extinction curves. With each extinction curve, attenuation curves 
were derived for three dust geometries (cloudy, dusty, and shell; see text). Bottom left: Models constructed with a power-law attenuation curve (r oc A" 7 ), with 
a fraction of starlight that is completely unobscured. Fractions of / = 0.012, / = 0.025, and / = 0.050 are shown. Bottom right: Power-law attenuation curves 
ranging from A~°- 4 to A~ 1,3 in steps of 0.3. 



cording to a two- comp onent dust model, the details of which 
are described in 33.21 The models differ principally in the 
attenuation curve associated with the diffuse cirrus dust. For 
each model, the sequence in colors is a sequence in increasing 
V-band optical depth, r^, associated with this diffuse compo- 
nent. The sequence increases from = 0.0 to = 1.5. 

In the upper left panel of Figure [3] we compare the sample 
of disk-dominated galaxies to MW attenuation curves with 
Ry = 2.0. The curves differ only in the treatment of the UV 
bump at 2 175 A. We consider curves w ith a UV bump char - 
acteristic of the average MW extinction (Car delli et al.ll 1989b . 
no UV bump, and a UV bump with strengt h equa l to 80% of 
the full value (B = 0.8 in the notation of 33.2. U . A model 
with a UV bump strength equal to 80% of the full value pro- 



vides an excellent fit to the data, while a model without a UV 
bump provides a poor match. For this model, the observations 
are reproduced for 0.2 < < 0.5. This range in tj is con- 
sistent with detailed radiative tran sfer models that consider 
inclination-dependent attenuation (Tuffs et aTll2004l) . 

While not shown, use of the Calzetti et al. attenuation curve 
for starbursts produces results qualitatively similar to the no 
bump curve, and therefore fails to reproduce the observed 
trends. Also not shown are models constructed with a MW 
curve with Ry = 3.1 (the average MW value). These mod- 
els also provide a poor match to the data, both for no UV 
bump and a bump with full strength {B = 1.0). A model with 
Ry = 3.1 and B = 0.6 provides a good fit to the UV data, al- 
though the best-fit values as inferred from the UV colors 
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cannot reproduce the observed optical colors (see $4.2i . For 
this reason, the Ry = 3.1 model is discarded. In contrast, the 
inclination-dependent UV and optical colors can be well-fit 
with the same values for both the UV and optical colors. 

We show model predictions for the attenuation curves of 
WGOO in the upper right panel of Figure [3] WGOO model 
predictions are presented for the 'homogeneous' local dust 
distribution, for all three large-scale geometrical configura- 
tions ('cloudy', 'dusty', and 'shell') utilizing both the MW 
and SMC extinction curves. For these models, unlike the oth- 
ers, the predictions are shown for optical depths ranging from 
0.0 to 50.0. The WGOO models with MW extinction all predict 
FUV-NUV colors too blue compared to the data. With SMC 
extinction the WGOO models qualitatively resemble the mod- 
els with the MW extinction curve without a UV bump shown 
in the upper left panel. The 'dusty' geometric configuration 
with SMC extinction is the model that begins to bend toward 
redder NUV-m colors. For this geometry the colors saturate 
because the attenuation curve becomes increasing greyer at 
high opacity (see Figure |2). One therefore cannot appeal to 
ever higher opacities to produce colors in better agreement 
with the data. The 'clumpy' local dust distribution predic- 
tions do not produce better fits to the data. It is of course pos- 
sible that a different geometrical configuration not considered 
by WGOO may yield better agreement with the data when the 
SMC extinction curve is used. Also, the spherical symmetry 
assumed in the WGOO models may be an important limitation 
for our purposes. It is however beyond the scope of this paper 
to consider such variations to the WGOO model. 

In the lower right panel we show results for power-law at- 
tenuation curves with an index that varies from 5 = -0.4 to 
5 = —1.3. These models all predict FUV-NUV colors far too 
red at a given NUV-m color. 

In the lower left panel we show results for a model con- 
structed that is identical to the power-law attenuation model 
with 6 = -0.7 except that a skin of completely unobscured 
stars is added with mass fraction of / = 1-5%. This model 
does somewhat better at describing the data than the pure 
power-law models in that the dependence of the FUV-NUV 
color on the NUV-m color becomes shallow at large opti- 
cal depths. The qualitative trends are a decent match to the 
data, but the quantitative predictions are much too blue in 
NUV-m. We will come back to this point below. Note also 
that the turnover in the model predictions at high opacities is 
generic to skin models. As the opacity increases, eventually 
the bulk population becomes heavily obscured and the (very 
blue) colors of the unobscured population dominates the in- 
tegrated color. If the opacity were increased beyond what is 
shown in the figure, the colors would eventually approach the 
Td = 0.0 colors. 

It is clear from Figure [3] that none of the commonly used 
attenuation curves provide a good description of the observed 
trends with inclination. This is so because the disk-dominated 
sample has FUV-NUV colors that vary only slightly with in- 
clination, while at the same time NUV-m varies by rs 0.5 mag. 
The implication is that the net attenuation in the FUV and 
NUV bands is comparable, i.e., Apuv ~ Anuv, and they both 
increase with increasing inclination. Most standard attenua- 
tion curves (e.g., a power-law, the form advocated by Calzetti 
et al. for starburst galaxies, a MW curve without a UV bump, 
and an SMC -based curve) do not satisfy this requirement. It is 
also difficult for geometrical effects to produce such a trend. 
The observed trend is however readily explained if the atten- 
uation curve contains a UV bump with strength equal to 80% 
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FIG. 4. — UV color-color diagram demonstrating the impact of various 
model assumptions. The default model (solid line), which assumes MW at- 
tenuation with a UV bump (B = 0.8) and Ry = 2.0 and a constant SFH, and 
the observational results (symbols), are repeated in all panels. Top left: Vari- 
ation in the V— band optical depth, t/x associated with young stars. Top right: 
Variation in the index, <5, of the attenuation curve for the young stars. Bottom 
left: Variation in the time during which young stars experience additional at- 
tenuation. This time is shown in the legend in units of log(i/yr). Bottom right: 
Variation in the SFH. Models are shown where a smaller fraction of mass is 
formed in a constant mode of SFH (c = 0.4 and c = 0.2); the remaining mass 
is assumed to have formed instantly at t = 0.0. Notice the different axis scales 
compared to Figuref5] 

of the canonical MW value, and with Ry = 2.0. In this case 
the net attenuation in the FUV and NUV bands is comparable, 
as evidenced by the top left panel of Figure [3] 

4. 1 . 1 . Assessing model assumptions 

In this section we consider the consequences of model vari- 
ations beyond the attenuation curve for diffuse dust. In Figure 
2] we explore the impact of several important model assump- 
tions on the UV colors. The default model here adopts the 
same underlying stellar population as before (constant SFH, 
solar metallicity), along with a MW attenuation curve with 
Rv = 2.0 and B = 0.8, tj, c = 0.5, a power-law attenuation with 
8 = —1.0 for the young stars, and a transition time between 
young and old stars of t = 10 7 yr. Variations to this default 
model that are shown in Figure [4] include the following: no 
additional attenuation around young stars (7j c = 0.0), twice 
the attenuation around young stars (jbc = 1 -0), an attenuation 
curve for young stars with an index of -0.7 and -1 .3, a transi- 
tion time between young and old stars, log(f esc /yr), of 6.7 and 
7.5, and SFHs in which only a fraction of stars are formed in a 
constant mode of star formation (c = 0.4 and c = 0.2), the rest 
having formed instantly at t = 0.0. 

In all cases the modifications simply result in a shift in the 
model locus in the FUV-NUV vs. NUV-m plane. There- 
fore, varying these model parameters allows us the freedom 
to move the model predictions in Figure[3]horizontally and/or 
vertically by a few tenths of a magnitude. Such shifts do not 
change our conclusions because the attenuation models that 



EVIDENCE FOR THE 2 175 A DUST FEATURE 



9 



-18.8 
-19.0 

-19.2 

-19.4 

-19.6 

-19.8 
-20.0 



• SDSS+GALEX data 
O model photometry 



face-on 
(0.9<b/a<1.0) 



edge-on 
(0.3<b/a<0.4) 



MW: R v =2.0, B=0.8 
i i i i i i i_ 



trv" 



♦ 



J*? 



i* 

II : 



MW:R V =3.1,B=0 
MW: R v =2.0, B=0 

j i i i i i i_ 



2000 4000 6000 
MA) 



8000 10000 2000 



4000 6000 
MA) 



8000 10000 



FIG. 5. — Average SEDs of face-on and edge-on star-forming galaxies with 9.5 < log(M/Mo) < 10.0 (solid symbols). Errors on the average SEDs are smaller 
than the symbol sizes. These data are compared to models with a constant SFH, Z = Zq, and a MW-like attenuation curve. Left panel: Our best-fit model 
attenuation curve with Rv = 2.0 and B = 0.8. The open symbols indicate the corresponding model magnitudes computed from the model spectra. Right panel: 
Model attenuation curves with R v = 2.0, B = 0.8 and R v = 3.1, B = 0.6. 



fail do so in their dependence of FUV-NUV on NUV-m. In 
other words, allowing the freedom to shift the models arbitrar- 
ily does not produce significantly better fits to the data, except 
for the skin models and the B = 1 .0 UV bump model; see be- 
low. Modifying the 'zero-point' of the models will however 
clearly change the quantitative conclusions drawn from this 
figure, such as the average SFH and total dust opacity, but 
that is to be expected when only considering UV colors. 

The only cases where a shift in the NUV-m color may re- 
sult in a better model fit is for the skin models and the B =1.0 
bump model. For the skin model, moving the model NUV-m 
colors redward by w 0.4 mag would yield an improved agree- 
ment between the skin models and the data. Such a shift can 
be induced by a SFH where only 20% of the mass forms in 
a constant mode (c = 0.2), the rest having formed instantly at 
t = 0.0. Below we show that such SFHs can be ruled out by the 
optical colors of our disk-dominated sample. For the B = 1 .0 
bump model, a shift redward in FUV-NUV of « 0.1 mag, in- 
duced for example by choosing f esc = 10 7 5 yr, would yield a 
good match with the data. We leave this as an open possibility, 
noting that such a shift would only strengthen our conclusion 
that a strong 2 175 A absorption feature was required to match 
the observations. 

4.2. Constraints from optical colors 

We now turn to optical colors in order to gain additional in- 
sight into the dust and stellar population content of our sam- 
ple. In Figure [5] we show the average spectral energy distri- 
butions (SEDs) for galaxies in our highest and lowest incli- 
nation bins. SEDs were constructed from GALEX and SDSS 
photometry. We also include in this figure our favored model, 
which has a constant SFH, Z = Z©, and a MW dust attenuation 
curve with R v = 2.0 and B = 0.8 (models with R v = 2.0, B = 0.0 
and Ry = 3.1,B = 0.6 are also shown for comparison). The 
models for the face-on and edge-on galaxies differ only in 
the V-band optical depth associated with diffuse dust, tj. 



Clearly our favored model matches the data well over the in- 
terval 1500A< A < 9000A, providing further support for the 
assumptions made in our favored model. The strong depres- 
sion in the SED of edge-on galaxies at A rs 2200A is striking, 
and strongly suggestive of the presence of the UV bump. 

The disagreement between the models and data in the 
z-band is noteworthy. Unfortunately, no significant conclu- 
sions can be drawn from this disagreement because the SED 
of star-forming galaxies at A > 7000A is sensitive to the 
prese nce of thermally -pulsating asymptotic giant branch stars 
(e.g.,|M a raston 2005). The model ing of this stellar evolution - 
ary phase is very uncertain (e.g.. IConroy et all 120091 l2010h . 
and minor modifications to the bolometric luminosity and 
temperature of these stars can readily account for the dis- 
agreement seen in the figure. 

In order to further assess our model assumptions, in Fig- 
ure|6]we consider optical color-color plots. In this figure the 
SDSS photometry is compared to models that differ in their 
attenuation curve associated with diffuse dust (left panel), and 
models that differ in their assumed SFH. Note that the model 
ugriz photometry is insensitive to the presence of a UV bump, 
and so we do not distinguish between MW attenuation with 
or without a UV bump in this figure. 

The variations in SFH result in significant differences in 
these ugr color-color plots. Notice that the bluest points in 
both the c = 0.4 and c = 0.2 SFH models are redder than the 
bluest data points. This implies that these SFHs cannot pos- 
sibly fit our disk-dominated sample (recall that a single SFH 
must be able to fit the colors for all inclination bins, since SFH 
will not vary with inclination). As mentioned in $2] this is not 
surprising in light of results from detailed modeling of the 
spectra of galaxi es with masses comparable to those consid - 
ered herein (e.g., ISahm et al.ll2007b iBrinchmann et al.ll2004h . 
We note in passing that the SFH models considered herein 
produce UV and optical colors that are almost indistinguish- 
able from r-model SFHs with r = 5 and 8 Gyr corresponding 
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FIG. 6. — Optical color-color plots as a function of inclination compared to 
various models. Left panel: Comparison of MW (with Ry = 2.0 and B = 0.8) 
and power-law (with <5 = —0.7) attenuation curves. Here the SFH is constant 
(c = 1.0). Each model is a sequence in optical depth associated with diffuse 
dust, in steps of A(r,j) = 0.3. Right panel: Models with the same MW atten- 
uation curve but varying SFH. The SFHs are varied such that a fraction, c, 
of the mass is formed in a continuous mode of star formation, with the rest 
having formed instantly at z = 0.0. 



to c = 0.2 and c = 4. 

As discussed in £ 14.1.11 the skin model shown in the lower 
left panel of Figure [3] could be brought into agreement with 
the data if a SFH with c = 0.2 were adopted. While not shown, 
the optical colors for the skin model are identical to those in 
Figure[6]for the MW curve because the flux contributed by the 
skin is negligible in the optical. The results from this figure 
therefore allow us to confidently rule out the skin models be- 
cause the SFH necessary to reconcile the model with the UV 
colors grossly fails to m atch the optical colors. 

As mentioned in 14.11 a model with Ry = 3.1 and B = 0.6 
also provided a good fit to the UV colors. However, as shown 
in Figure [5] this model cannot simultaneously match the en- 
tire UV through near-IR SEDs of star-forming galaxies. The 
disagreement is particularly striking for the NUV-m colors, 
where the constraints on low values for Ry are strongest. 

4.3. Trends with stellar mass and metallicity, and variation 
in data analysis 

We now return to the UV colors. In Figure [7] we explore 
the sensitivity of the results shown in Figure [3] to the sam- 
ple selection and method of data analysis. In this figure we 
consider galaxies both more and less massive than our fidu- 
cial sample, and galaxies split according to their gas-phase 
metallicity. In the latter case we consider a wider bin in 
mass, 9.0 < log(M/M ) < 10.0, in order to increase statis- 
tics. In the panel showing the dependence on metallicity, we 
also show the change in model colors accompanying a change 
in metallicity of the underlying stellar population from 0.4Z Q 
to Z©. Model predictions for the default set of assumptions 
are included in this figure, where we have adopted a MW at- 
tenuation curve with Ry = 2.0 and with UV bump strengths 
B = 0.0,0.8,1.0. These model predictions are identical to 
those shown in the upper left panel of Figure [3] 

The observed trends are generally as expected: more mas- 
sive and more metal-rich galaxies are redder in both NUV-m 
and FUV-NUV, although the differences are not large. There 
is tentative evidence that the more massive galaxies deviate 
from our favored model in the sense that their FUV-NUV col- 
ors are w 0.1 mag bluer than expected for edge-on galaxies. 
While we are hesitant to interpret this trend, if real it could 
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FIG. 7. — Color-color diagram comparing a variety of different selection 
criteria and methods of data analysis. Models with a MW attenuation curve 
with varying UV bump strength (B = 0.0.0.8, 1.0) are included for compari- 
son (solid lines; these models are identical to the models shown in the upper 
left panel of Figure f3j. In each panel the fiducial sample is shown as dia- 
monds. Notice the different axis scales compared to Figure [5] Upper left: 
Samples selected with stellar masses 9.0 < log(A/ /Mq) < 9.5 and 10.0 < 
log(M/M Q ) < 10.5. Upper right: Samples selected with 12+[0/H]< 9.0 and 
12+[0/H]> 9.0 drawn from a parent sample with 9.0 < log(M/M Q ) < 10.0. 
The arrow indicates the change in colors accompanying a change in Z of 
the underlying stellar population from 0.4Zq to Zq. Lower left: FUV-NUV 
color computed after removing the 5% bluest pixels. Lower right: Compari- 
son between colors of the fiducial sample computed by averaging magnitudes 
and averaging fluxes. 



indicate a stronger UV bump in more massive galaxies. 

The observed weak dependence on metallicity is generally 
consistent with being due to the underlying stellar population, 
under the assumption that variation in gas-phase metallicity is 
accompanied by a similar change in the stellar metallicity. 

In Figure|7]we also show the observed trend where we have 
computed the FUV-NUV color after masking the bluest 5% of 
pixels. This masking was done to test the extent to which the 
observed trend (namely, the weak trend of FUV-NUV vs. in- 
clination) is due to a few very blue and concentrated regions 
of star formation. The NUV-m color was not re-computed 
with this masking because of the complications in translating 
the masks in the GALEX pixels into the SDSS pixels. The 
resulting FUV-NUV colors are on average » 0.1 mag redder 
than the colors computed from all of the pixels. The results 
from this exercise demonstrate that the observed very blue 
FUV-NUV colors even for the most inclined galaxies cannot 
be attributed to a few very blue concentrated star-forming re- 
gions. 

Finally, in this same figure we show how the average colors 
of our fiducial sample depends on the method of averaging. 
We show average colors computed by averaging fluxes (our 
standard a ppro ach) and by averaging magnitudes. As men- 
tioned in q2.ll in the absence of dust, averaging fluxes is the 
correct approach. In the presence of dust, there is no clearly 
preferred method, although we note that the attenuation by 
dust is averaged in a more meaningful way when averaging 
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FI G. 8. — Total FUV attenuation as a function of FUV-NUV color (ana log of the IRX-/3 relation) . Lines show the relation for starburst galaxies (dashed 
/me; Kong et al. 2004). for galaxies in the SINGS survey (dot-dashed tine; Munoz-Mateos et al. 2009), and our best-fit model that contains a strong UV bump 
(solid line). For the starburst and SINGS samples, Apuv is estimated from Lir/Lfuv according to the relation provided by Buat et al. 12005). For our samples 
of disk-dominated galaxies (diamond, star, and square symbol s), Apuy is estimated as described in the text. We show results for three bins in stellar mass, 
as indicated in the legend. Small symbols are from the work of Johnson et al. 1 2007a, J07); larger connected symbols are median FUV-NUV colors from J07 
computed in bins of Apuv ■ 



magnitudes. Regardless, the resulting average inclination- 
dependent colors are largely unaffected by our choice of av- 
eraging. 

4.4. The IRX-P relation 

Recall that the IRX-/? relation compares the UV spectral 
slope, (3, to the ratio of total infrared to UV luminosity, the 
latter being a fairly direct probe of the total UV attenuation. 
In theory, an lRX-(3 relation therefore provides a simple es- 
timate of the net attenuation in the UV based solely on the 
UV slope (i.e., the FUV-NUV color). In practice, this proce- 
dure is complicated not only by the large observed scatter in 
the relation, but also by the fact that starburst and 'normal' 
star-forming galaxies follow different relati ons dMeurer et alj 
1 1999t iDale etaLll200 7 : Boissi er et al1l2007l). and the fact that 
the re lation depends on SFHfKong et al.l l2004t[Cortese etal] 
2008). 

In Figure [8] we translate our results into an analog of the 
IRX-/? relation. Specifically, we estimate the attenuation in 
the FUV, Apuv, of the observed galaxies from their g — r color 
by simply comparing the model colors to the observed col- 
ors in order to estimate the total net attenuation [i.e., Apuv oc 



ln(^ FUV /,F FUV ) where i and o denote the dust-free and dust 
reddened model spectra, respectively]. We obtain similar re- 
sults when using other colors to estimate Apuv- We then shift 
our estimated Apuv values down by 0.8 mag in order to match 
the locus of previous results (see below). The particular val- 
ues of Apuv will clearly depend on our model assumptions 
(e.g., the assumed SFH), but the general behavior of our re- 
sults in the IRX-/3 should be robust. Results are shown for 
three stellar mass bins. 

In this figure we include the IRX-/3 relation estimated 
for a sample of local starburst galaxies (Meurer et al.lll999t 
iKong e t al. 200 4J), and normal star-forming galaxies from the 
SINGS survey (iMunoz-Mateos et al.ll2009h . The TRX' por- 
tion of these relatio ns has been co nverted into Apuv via the 
formula provided in Buat et al] (120051) . 

The most striking conclusion to be drawn from Figure [8] is 
that our sample spans a large range in Apuv but a compara- 
tively narrow range in FUV-NUV color, i.e., our results imply 
an extremely steep IRX-/3 relation. The relation for our sam- 
ple of disk-dominated galaxies is much steeper than either 
the starburst or normal star-forming galaxy relation. This re- 
sult is driven entirely by the fact that while FUV-NUV varies 
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little in our sample, other colors, such as NUV-m, u — r, and 
g—r vary significantly. The trend for the highest stellar mass 
bin is notable in that it ap pear s to be steeper than our favored 
model. As mentioned in 8 14.31 if real this trend could point to 
a stronger UV bump in more massive galaxies. 

The IRX-/3 relation from the SINGS sample is rep- 
resentative of the entire SINGS sample, excluding only 
three starburst galaxies an d galaxies with FUV-NUV>0.9 
(Munoz-M ateos et al.ll2.009t) . It is therefore not surprising that 
the SINGS relation is much shallower than ours, because their 
sample contains g alaxies with a wide range in SFHs. Indeed, 
from Figure 3 in Munoz-Mateos et al. (2009) it is apparent 
that their TRX-/3 relation is driven largely by variations in 
SFH, in the sense that galaxies with redder FUV-NUV col- 
ors have formed the bulk of their stars at earlier epochs. In 
contrast, the IRX-/3 relation derived herein is characteristic 
of a single SFH (for each stellar mass bin), and therefore our 
relation is truly a sequence in dust opacity. 

In th is figure we also compare our results to lJohnson et alj 
d2007al J07) who computed the IRX-/3 relation for wlOOO 
galaxies using observations from GALEX and Spitze^ We 
show their results for all galaxies with D„4000 < 1.5, where 
D„4000 is a measure of the strength of the 4000A break. 
D„4000 is a measure of the age of the stellar population that 
is relatively insensitive to dust. This cut on D„4000 selects 
actively star-forming galaxies, similar to those in our own 
sample. For this sample we also show the median FUV-NUV 
color in bins of Apuv, where the FUV attenuation has be en 
estimated from IRX using the formula in Bu at et ail (|2005). 

The average relation from J07 is in excellent agreement 
with our observational results and also with our favored model 
that includes a strong UV bump in the diffuse dust attenuation 
curve, after we shift our results downward in Figure [8]by 0.8 
mag. We stress that the solid line in Figure[H]is not a fit to the 
average relation from J07. Rather, the agreement is remark- 
able independent confirmation of our best-fit model. We will 
discuss the implications of this result in the following section. 

The required shift of our results 0.8 mag downward to 
match the results from J07 deserves further comment. Note 
first that the results from J07 and those derived herein are not 
directly comparable, as J07 measure IRX while we only mea- 
sure UV through near-IR SEDs. Any attempted comparison 
between our work and J07 will therefore be model-dependent. 
A detailed analysis of this issue is beyond the scope of the 
present article. However, irrespective of our results, notice 
that the implied Apuv vs. FUV-NUV relation in the J07 sam- 
ple extrapolates to Apuv = 0.0 at a FUV-NUV color signifi- 
cantly redder than canonical dust-free solar metallicity stellar 
population predictions (FUV-NUV« 0). This indicates that 
there may be some modeling error in translating IRX (which 
J07 actually measure ) into Apuv (wh ich was derived from IRX 
via the relation from Buat et al. 2005), at least at low attenua- 
tion. 

5. DISCUSSION 

5.1. An emerging physical picture, and implications 

In the previous section we demonstrated that the 
inclination-dependent UV colors of a mass-selected sample 
of disk-dominated galaxies are best explained by a MW at- 

6 The data from J07 shown in Figure[8]are revised quantities based on the 
same sample as in J07. The most significant change with respect to the data 
presented in J07 is the use of updated GALEX photometry, which results in 
bluer FUV-NUV colors by » 0.05 - 0. 1 mag. 
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FIG. 9. — Attenuation curves normalized to the V— band. The dotted curves 
show MW attenuation curves with Ry = 2.0, both for a standard UV bump 
strength and no UV bump (B = 1 .0 and B = 0.0, respectively). The thick solid 
line is our best-fit attenuation curve with Ry = 2.0 and a UV bump strength 
equal to 80% of the standard MW value. The thin solid line is the standard 
MW curve (i.e., with Ry = 3.1 and B = 1.0). The points show the actual 
constraints on the attenuation curve provided by the broadband photometry. 
Transmission curves with arbitrary normalization are included for the FUV, 
NUV, u, and g filters. 



tenuation curve with Ry = 2.0 that includes a prominent UV 
bump at 2 175 A. The strength of this bump is approximately 
80% as strong as the average UV bump observed in the MW. 
This result is summarized in Figure [9] In this figure we show 
our favored attenuation curve and compare this to the standard 
MW curve (Ry = 3.1). We also show the actual constraints on 
the attenuation provided by the FUV, NUV, u, and g filters. 
It is clear from this figure that while our results are consis- 
tent with a strong UV bump, we cannot claim that our results 
require a UV bump, since we rely solely on broadband pho- 
tometry over a narrow redshift range. However, inspection of 
Figure [5] strongly suggests that the UV bump is responsible, 
as we know of no other model variation capable of reproduc- 
ing the depression in the average SED of edge-on galaxies at 
A « 2200A. 

Our results imply that previous work attempting to mea- 
sure the SFR in low redshift galaxies has underestimated the 
amount of attenuation in the near-UV. From inspection of Fig- 
ure [9] we can see that typical attenuation in the NUV band is 
underestimated by approximately the value of the V-band op- 
tical depth, Ty, when comparing standard attenuation curves 
(power-law, starburst, and average MW) to our favored model. 
For moderately dusty galaxies the underestimation may there- 
fore be substantial. The bias is such that SFRs based primarily 
on the near-UV will be lower than the intrinsic SFR. Addi- 
tional work will be required to understand in detail how the 
effects of a strong UV bump propagate into the derived phys- 
ical properties of galaxies. 

From this figure we can also understand qualitatively why 
our disk-dominated sample favors an attenuation curve with 
Rv = 2.0, as opposed to the canonical value of Ry = 3.1 pre- 
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ferred in the MWQ Decreasing Ry from 3.1 to 2.0 results in 
substantially more attenuation between the u and NUV filters, 
and therefore a much redder NUV-m color. The relative at- 
tenuation between FUV and NUV, and also between u and 
g is does not change substantially between the R v = 3.1 and 
Rv = 2.0 curves, and so these colors change by much less than 
NUV-w. Therefore, it is principally the constraint from the 
NUV-m color that drives the requirement for Ry ~ 2.0. 

If our interpretation of the observed trends is correct, this 
would constitute the first detection of the UV bump in the 
attenuation properties of galaxies at low redshift. The only 
other detection of the UV bump in the att enuation curves 
of galaxies was reported by iNoll e t al. (2009), who consid- 
er ed restframe UV spectra of star-forming galaxies at z ~ 2. 
In IConrovl ( 120101) the observed B-R colors of star-forming 
galaxies at 0.7 < z < 1.4 were used to constrain the average 
strength of the UV bump. A bump as strong as observed in 
the MW was ruled out for Ry = 3.1, as we find herein, but 
a UV bump strength needed to match our sample (B w 0.8 
with Ry = 2.0, which, in the UV, is comparable to B pa 0.6 
with Ry = 3.1) could not be ruled out. More detailed model- 
ing of the z ~ 1 population will be required to place stronger 
constraints on the presence of the UV bump. Luckily, the nec- 
essary data already exists. 

It is important to understand why evidence of the UV bump 
in the attenuation curve of galaxies has proved elusive until 
now. Many authors have previously attempted to model the 
UV colors of galaxies with a standard MW attenuation curve 
(i.e., with Ry = 3.1 and B = 1.0). In all cases the MW curve 
was found to be a poor fit, as we find herein. The reduction in 
the strength of the UV bump to B = 0.8, accompanied by our 
choice of Ry =2.0, is essential for our attenuation models to fit 
the data. Recall again that an attenuation curve with Ry = 2.0 
and B = 0.8 yields FUV-NUV colors comparable to an atten- 
uation curve with Ry = 3.1 and B = 0.6. In other words, if 
considering only the UV and Ry = 3.1 MW-like attenuation 
curves, a significant depression of the UV bump is required. 
Clearly, ruling out the standard MW curve does not imply that 
observed galaxies show no evidence for a UV bump. In ad- 
dition, considering the UV colors of galaxies as a function of 
inclination allowed us to confidently separate stellar popula- 
tion effects from dust attenuation effects. 

These results must somehow be reconciled with the well- 
known fact that lo cal starburst galaxies show no evidence for 
a UV bump ("e.g.. ICalzetti et al]|1994l) . It has been conjec- 
tured that either metallicity, intensity of the interstellar radi- 
ation field, or the star-dust geomet ry modulates the s trength 
of the U V bump. The dust mod el of Sil va et aU (119981) . as ex- 
plored in Granato et al. ( 2000), provides a plausible path to- 
ward reconciling our results with those for the starburst galax- 
ies. Granato et al. showed that an underlying MW extinction 
curve (applied to dust in the diffuse ISM) can give rise to a 
starburst attenuation curve similar to that observed by Calzetti 
et al., and can also give rise to an attenuation curve with a sig- 
nificant UV bump for 'normal' star-forming galaxies. In their 
model the difference is due to geometry: for the starbursts the 
UV dust attenuation is dominated by dust within the molec- 
ular clouds in which the young stars are embedded, because 
the fraction of young stars is high in starbursts. The molec- 

7 Of course, as emphasized throughout, attenuation and extinction are con- 
ceptually different, and so there is no a priori reason for the attenuation curve 
that best describes our sample to be similar to the average extinction curve 
measured for the MW. 



ular clouds are optically thick, and so the resulting attenua- 
tion curve is governed primarily by the wavelength-dependent 
fraction of the intrinsic flux emitted by young stars. In the 
case of normal star-forming galaxies, a significant fraction 
of UV photons are emitted by stars that have left their birth 
clouds, and so the only attenuation they suffer is due to the 
optically thin diffus e ISM, which may imprint a UV bump. 

More recently, dPanuzzo et alj|2007h have used the Silva et 
al. dust models to interpret the IRX-/3 relation for normal 
SF galaxies and conclude that the data are consistent with a 
UV bump in the extinction curve. These authors highlight the 
important role of a ge-dependent exti nction on the resulting 
attenuation curve. Ilnoue et all d2006l) consider the effects of 
dust scattering in detail and co nclude that the data are con- 
sistent with a UV bump. Finally, Burgarella et al. (2005|) take 
an empirical approach in modeling UV attenuation similar to 
what we do herein, and conclude that the observed IRX-/3 re- 
lation for normal SF galaxies is consistent with a UV bump 
with half the strength of the UV bump in the MW extinction 
curve. All of these studies provide additional support to our 
principle conclusion regarding the presence of a UV bump in 
the attenuation curve of normal SF galaxies. 

While the Silva et al. model is capable of explaining the 
absence of a UV bump in starburst galaxies and its presence 
in normal star-forming galaxies by geometrical effects, other 
possibilities remain. The most plausible alternative explana- 
tion is the preferential destruction of the grains responsible for 
the UV bump in the harsh interstellar radiation characteristic 
of starburst galaxies. PAH emission in the IR is known to be 
modulated by the strength of the radiation field in the sense 
that harder radiation fields result i n dim i nished (or even neg- 
ligible) PAH emi ssion (e.g.. IVoiil[l992l: iMadden et al J 120061: 
Smith et al. 2007). These properties of PAHs can provide a 
natural explanation for the variation in UV bump strength, if 
PAHs are responsible for the absorption at 2 175 A. More di- 
rectly, the UV bump is absent in the extinction curve probed 
by four out of five sightlines in the SMC dPeil 1 19921). which 
may be due to the harsher radiation field there (Gord onet alJ 
2003). Metallicity is an unlikely explanation for the differ- 
ence between starburst and normal galaxies because the gas- 
phase metallicities of these two types of galaxies are not sub - 
stantially different ( ICalzetti et al.|[l994l:lTremonti et alj2004t) . 
Now that we have identified samples of galaxies that plausibly 
contain a UV bump, each of these scenarios can be tested. 

In Figure[8] we showed that the average relation between to- 
tal UV attenuation and UV spectral slope is extremely steep in 
the sense that a narrow range in observed spectral slope (i.e., 
FUV-NUV color) corresponds to a large range in UV attenua- 
tion. These results extend and strengthen previous results that 
also fo und very steep relations (e.g., Bell 2002; Johns on et alJ 
2007b). In our interpretation, the steepness of this relation is 
due to the properties of the average attenuation curve for a 
fixed average SFH. In essence, the attenuation curve flattens 
considerably over the wavelength range probed by the FUV 
and NUV filters (1500A< A < 2300 A), and so large variation 
in the total attenuation produces little change in FUV-NUV 
colors. Appeal to the UV bump at 2175A provides a nat- 
ural explanation for the observed flattening in the attenuation 
curve. The shallower IRX-/3 relations found in previous work 
arises from the fact that previous relations were constructed 
from samples of galaxies spanning a wide range of SFHs. 

We have presented average UV colors as a function of incli- 
nation and UV attenuation. For individual galaxies, one can 
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FIG. 10. — Restframe UV spectra of two models that differ only in the 
strength of the UV bump at 2 175 A. One model has no UV bump while the 
other has a UV bump with strength 80% of the average MW strength. Each 
model is a sequence in optical depth characterizing the diffuse dust compo- 
nent: Trf = 0.0,0.3,0.9. The spectra are normalized to the flux at 3000A. 

expect variations about these average relations due to vari- 
ations in either the star-dust geometry, SFH, metallicity, or 
dust-to-gas ratio. Indeed, within a given inclination bin the 
spread in FUV-NUV colors is large (cr « 0.5 mag). We spec- 
ulate that variation in the SFH over the last ~ 10 8 yr and vari- 
ation in the star-dust geometry are the two largest sources of 
scatter. In addition, we clearly detect systematic changes in 
the FUV-NUV color with stellar mass, which may either be 
due to a varying attenuation curve or possibly a varying aver- 
age SFH with galaxy mass. 

In light of the steepness of the TRX-/3 relation for a given 
SFH as measured herein, its dependence on the SFH, and the 
expected galaxy-to-galaxy scatter due to variations in other 
physical properties, we conclude that the uncertainties in es- 
timating UV attenuation based on the observed UV slope are 
very large and believe that this calls the utility of the IRX-/? 
relation into question. This conclusion applies when j3 — the 
UV slope — is measured over a wavelength interval includ- 
ing the UV bump. If this spectral range is avoided then the 
IRX-/3 relation may yield more reliable results. Care must 
clearly be taken when estimating j3 from either spectroscopy 
or photometry. 

5.2. Future directions 

The results contained in this work present clear directions 
for future work. First, moderate resolution spectra in the UV 
would provide conclusive evidence for the UV bump. In Fig- 
ure [10] we show the restframe UV spectra for our default stel- 
lar population model (constant SFH and solar metallicity), for 
two dust attenuation curves that differ only in the strength of 
the UV bump. We show models both for a bump strength of 
80% (the favored dust model) and zero bump strength. Spec- 
tra are shown for three values of the V-band optical depth 
characterizing the diffuse dust: Td = 0.0,0.3,0.9. The treat- 
ment of dust surrounding young stars is constant and the same 
as previous sections. For reference, the average inclination- 



dependent colors of our disk-dominated sample can be repro- 
duced in our favored model for 0.2 < Td < 0.5. 

This figure demonstrates that UV spectra of moderately 
dusty galaxies would provide a clear test of our preferred 
model. Our sample of highly inclined disk-dominated star- 
forming galaxies would be ideal for these purposes. Narrow 
band photometry around 2 175 A could also potentially reveal 
the unique signature of the UV bump. The required filters are 
already on-board the Hubble Space Telescope. 

If the UV bump is indeed responsible for the observed 
trends, the next task will be to identify whether or not PAHs 
are responsible for the absorption. Recall that the identity of 
the grain(s) responsible for the UV bump is still a mystery, 
and, while PAHs are a plausible candid ate, other possibilitie s 
remain (including, e.g., graphite; Draine & Ma lhotralll993l) . 
If PAHs are responsible, then one would expect to uncover 
a strong correlation between PAH emission in the infrared 
(traced for example by the 8/xot chann el on the Spitzer Space 
Telescope; Munoz-Mateos et al. 2009) and a strong UV bump 
(traced for example by the FUV-NUV color). 

An especially promising direction would be to consider 
galaxies with gas-phase metallicities both above and below 
12+[Q/H|= 8 1. B ased on data from the SINGS survey, 
Draine ~etaT] d2007l) found that galaxies above this metallicity 
contained a wide range of PAH emission strengths, whereas 
below this value galaxies were deficient in PAH emission. We 
would therefore expect the trends of FUV-NUV with inclina- 
tion to change qualitatively below this metallicity scale, again 
as suming that PAHs are t he cause of the UV bump. According 
to iTremonti et alj (|2004), galaxies with such low gas-phase 
metallicities have stellar masses log(M/M Q ) < 8.5. Irregular 
galaxies are more common than disk-dominated galaxies at 
this mass scale, and so it may be difficult to separate trends 
with geometry from trends with PAH abundance across this 
mass scale. Nonetheless, with detailed modeling the PAH hy- 
pothesis may be tested with the approach outlined above. 

6. SUMMARY 

In this work, we measure the UV and optical colors of disk- 
dominated galaxies at 0.01 < z < 0.05 as a function of their 
inclination. These inclination-dependent colors provide pow- 
erful constraints on dust attenuation models because the ob- 
served trends cannot be caused by variations in the underlying 
stellar populations. The observed inclination-dependent FUV- 
NUV and NUV-m colors cannot be explained by common 
dust attenuation curves, including power-law models, stan- 
dard MW and SMC cur ves, nor the att enuation curve mea- 
sured for starbursts dCalzetti et al.ll 19941) . After considering a 
number of model variations, we are led to conclude that the 
strong dust feature at 2 175 A (the UV bump) seen in the MW 
extinction curve is responsible for the observed trends. In or- 
der to simultaneously fit the UV and optical colors, we require 
an attenuation curve analogous to the MW extinction curve 
with R v = A v /E(B-V) « 2.0 and a UV bump strength equal 
to 80% of the standard MW value. Previous work failed to 
recognize the importance of the UV bump because only stan- 
dard MW curves were considered (Ry = 3.1, full UV bump 
strength), and these consistently fail to match observed trends. 

We also construct an analog of the TRX-/3 relation and find 
that the inferred relation between UV slope and FUV atten- 
uation is extremely steep. Our derived relation agrees well 
with recent work that combin es UV and IR phot ometry for 
«1000 low redshift galaxies dJohnson et al. 2007^), but is 
much steeper than the starburst relation. Our results provide 
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a natural interpretation of the IRX-/3 relations found in pre- 
vious work. Namely, for a given SFH, the IRX-/3 relation 
will be very steep because of the properties of the attenua- 
tion curve. Varying the SFH will produce primarily horizontal 
shifts in this relation. When considering heterogeneous sam- 
ples of galaxies with a range of SFHs, the IRX-/3 relation will 
on average appear shallow, due to the range in SFHs, but will 
show tremendous scatter, due to the attenuation curve. From 
the steepness of this relation for a given SFH, we conclude 
that it is not possible to accurately correct for dust attenuation 
— even in an average sense — based solely on the UV slope. 

Confirmation of our interpretation must await UV spectro- 
scopic follow-up. If confirmed, these results will pave the way 
toward finally identifying the grain population responsible for 
the 2 175 A extinction feature. In addition, these results neces- 
sitate a more nuanced interpretation of the UV spectral slopes 
of star-forming galaxies. 
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APPENDIX 

THE MILKY WAY EXTINCTION CURVE FOR ARBITRARY UV BUMP STRENGTH 

In this Appendix, we provide formulae to construct an extincti on curve analogous to the average MW extinction curve, but for 
arbitrary UV bump strength. We adopt the ICardelli et alj ([1989) parameterization of the MW extinction and only modify their 
formulae where necessary. Below, B will be the single parameter characterizing the strength of the UV bump at 2 175 A. A value 
of B = 0.0 results in no UV bump, while a value of B = 1.0 returns the standard MW UV bump strength. An IDL routine that 
computes the extinction curve for arbitrary Ry and B is available upon request from the authors. 

Cardelli et al. parameterize the MW extinction curve in the following way: 



A(X)/A v = a(x) + b(x)/R, 



(Al) 



where x has units of [im ,R = E(B-V)/Av, A(X) is the wavelength-dependent extinction in magnitudes, Ay is the extinction in 
the V-band, and a(x) and b(x) are specified in the following piece-wise continuous way: 
In the infrared: 0.3/im" 1 < x < l.l/ira -1 , 



a(x) = 0.574x L61 , 
b(x) = -0.527x lbl . 
In the optical/near-IR: 1 . 1 /im -1 < x < 3 .3/xm and y = x - 1 .82, 

aW=l+0.177y-0.504y 2 -0.0243y 3 + 0.72l/+0.0198y 5 -0.775/ + 0.330y 7 , 
Z?(x)=1.413y + 2.283y 2 +1.072y 3 -5.384y 4 -0.622y 5 + 5.303y 6 -2.090y 7 . 
In the near/mid-UV: 3.3/im" 1 < x < 5.9fim~ l , 

f a = f—\ (-0.0370 + 0.0469fi-0.601B/fl + 0.542/7?), 



a(x)= 1.752 -0.3 16x- 



6O) = -3.09+1.825x + 



0.104B 



(x-4.67) 2 + 0.341 

1.206B 
(x- 4.62) 2 + 0.263' 



+ fa, 



In the far-UV: 5.9^ < x < 8.0/im" 1 , 



f a = -0.0447 (x- 5 .9) 2 - 0.00978 (x-5.9) 3 , 
f b = 0.213 (x-5.9) 2 +0.121 (x-5.9) 3 , 



(A2) 
(A3) 

(A4) 
(A5) 



(A6) 
(A7) 
(A8) 



(A9) 
(A10) 
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0.104B 

a (x)=1.752-0.316x- (x _ 467)2 + () 34i +/a , (All) 

^) = -3.09 + 1.825x + (x _ 4 1 6 g 2 6 ; a263 +/ , (A12) 



